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Abstract 

We study multivariate trigonometric polynomials, satisfying a set 
of constraints close to the known Strung-Fix conditions. Based on the 
polyphase representation of these polynomials relative to a general 
dilation matrix, we develop a simple constructive method for a spe- 
cial type of decomposition of such polynomials. These decompositions 
are of interest to the analysis of convergence and smoothness of mul- 
tivariate subdivision schemes associated with general dilation matri- 
ces. We apply these decompositions, by verifying sufficient conditions 
for the convergence and smoothness of multivariate scalar subdivi- 
sion schemes, proved here. For the convergence analysis our sufficient 
conditions apply to arbitrary dilation matrices, while the previously 
known necessary and sufficient conditions are relevant only in case of 
dilation matrices with a self similar tiling. For the analysis of smooth- 
ness, we state and prove two theorems on multivariate matrix sub- 
division schemes, which lead to sufficient conditions for C 1 limits of 
scalar multivariate subdivision schemes associated with isotropic di- 
lation matrices. Although similar results are stated in the literature, 
we give here detailed proofs of the results, which we could not find 
elsewhere. 

MSC 41A30, 42C40 



1. Introduction 

In this paper we study multivariate trigonometric polynomials (masks), satis- 
fying a set of constraints, which is close to the known Strung-Fix conditions. 

1 This research was supported by The Hermann Minkowski Center for Geometry at Tel- 
Aviv University. The second author is also supported by DFG Project 436 RUS 113/951 
and Grant 09-01-00162 of RFBR. 
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We are interested in decomposing such a polynomial multiplied by a factor 
of the form 1 — e 2m ( x < e j)^ into a sum of d trigonometric polynomials, satisfy- 
ing the same type of constraints, but of one order less, each multiplied by a 
different factor of the form 1 — e 2m ( M * x > e k) _ 

This type of decompositions is of interest in the analysis of convergence 
and smoothness of multivariate subdivision schemes. The algebraic approach 
to the analysis of such schemes is investigated in [1] relative to the dilation 
matrix 21, and in [4], [8], [9], relative to general dilation matrices. 

The analogous decompositions for multivariate algebraic polynomials are 
based on the rich structure of ideals of multivariate polynomials. The trigono- 
metric decompositions are based on the rather simple idea of the representa- 
tion of a trigonometric polynomial in terms of its unique polyphase trigono- 
metric polynomials relative to a dilation matrix [10]. 

The outline of the paper is as follows: in Section 2 we introduce our 
notation and the set of constraints. We also bring basic results related 
to polyphase representations of trigonometric polynomials. Properties of 
trigonometric polynomials satisfying the set of constraints of order zero rel- 
ative to a general dilation matrix are derived in Section 3, in terms of their 
polyphase trigonometric polynomials. 

These results are used in Section 4 to prove the decomposition of trigono- 
metric polynomials, satisfying the set of constraints of order zero. An algo- 
rithm for the computation of this decomposition is presented. The decom- 
position for order zero is used to derive that for order one. The general case 
of order greater than one, which is much more involved, is also investigated, 
and an algorithm for this case is given. 

Section 5 gives a characterization of a trigonometric polynomial satisfying 
a set of constraints of a general order, in terms of the values at the origin 
of the derivatives of its polyphase trigonometric polynomials. This result 
provides a tool for the construction of such polynomials. 

Applications of the decompositions of Section 4 to the convergence and 
smoothness analysis of multivariate subdivision schemes associated with gen- 
eral dilation matrices are presented in Section 6. The application of the 
decompositions to the convergence analysis is based on a new result giving 
sufficient conditions for the convergence of a scalar multivariate subdivision 
scheme associated with an arbitrary dilation matrix. The result in the liter- 
ature giving necessary and sufficient conditions for convergence, is limited to 
dilation matrices having a self similar tiling [4]. For the analysis of smooth- 
ness, we state and prove two theorems on multivariate matrix subdivision 
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schemes. First we give a detailed proof of a sufficient condition for C 1 limit 
functions of such a scheme, associated with an isotropic dilation matrix. Al- 
though this condition is stated in [9], the sketch of a proof given there, is 
valid only for dilation matrices which are multiples of the identity matrix. 
Then we state and prove a condition, in terms of the decompositions, that 
allows us to check the sufficient condition for C 1 . 



2. Notation and preliminary information 

Let N be the set of positive integers, W d denotes the ci-dimensional Euclidean 
space, x = (xi, . . . , xj), y = (yi, ■ ■ ■ ,Ud) are its elements (vectors), (x, y) = 
XiVi + ■ ■ - + %dyd, M = m &x \xj\, e,- = (0, . . . , 1, . . . , 0) is the j-th unit vector 

j=l,...,d 

in R d , = (0, . . . , 0) G R d ; Z d is the integer lattice in R d . For x, y G R d , we 
write x > y if xj > yj, j = 1, . . . , d; Z d + = {x G Z d : x > 0}. If a, (3 G Z^, 

d / \ . d d 

a,b G M. d , we set a\ = ]J aijl, (J) = p^zm. , « & = II H = E a j, 

j=l V / j=\ j=\ 

D a f = — Q adx ; 5 a i, denotes the Kronecker delta; z := (z\, . . . , z d ), where 

Zk ;= e 2^lr = fN ;= {/ = {/Q}a6zd . fa e ^^f^ < 

WfWeg = ll/lloo := sup aGZd \f a \, £oo := 4,- 

Let M be a non-degenerate d x d integer matrix whose eigenvalues are 
bigger than 1 in module, M* is the transpose of M, Id and denote re- 
spectively the unit and the zero d x d matrices. We say that the numbers 
k, n G Z d are congruent modulo M (write k = n (mod M)) if k — n = M£, 
£ G Z d . The integer lattice Z d is splitted into cosets with respect to the intro- 
duced relation of congruence. The number of cosets is equal to | det M\ (see, 
e.g., [f f, p. 107]). Let us take an arbitrary representative from each coset, 
call them digits and denote the set of digits by D(M). Throughout the paper 
we consider that such a matrix M (dilation matrix) is fixed, m = \ detM|, 
D(M) = {s , . . . , Sm _x}, D(M*) = {s* , . . . , s = s* = 0, r k = M~ l s k , 

k — 0, . . . , m — 1. 

If A is a N x N' matrix with entries a^, then denotes its n-ih 
Kronecker power, i.e., 

/ a n AW . . . a 1N .AW \ 

A^ : = 1, A® := A, A [n+1] : = 

\a N1 AW ... a NN ,AWj 
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Proposition A The matrix < -4= e 2nt( - rk > s i) } is unitary. In particular, 

I vm J fc,Z=0 

1 m— 1 

^Ve 2 *^)^ ffl . (1) 



fc=0 



A proof of this statement can be found in [6] or in [2]. 

We will consider 1-periodic trigonometric polynomials in d variables 

t(x) = ^t{n)e 2 ^ x \ 



where the set {n G Z d : t(n) ^ 0} is finite. 

For any trigonometric polynomial t, there exists a unique set of trigono- 
metric polynomials r u , v = 0, . . . , m — 1, (polyphase functions of t) (see, 
e.g., [6]) such that 

m— 1 

= ^e^^V^M*^. (2) 

i/=0 

It is clear that 

Tu {x) = t( Mn + s v )e 2m{n > x) . 
nei d 

For any n G Z+, denote by the set of trigonometric polynomials t 
satisfying the following "zero- condition" of order n: D /3 t(M*~ 1 x)\ x=s = for 
all s G D(M*), s ^ 0, and for all f3 G Z%, [(3} < n. It is well known that 
this "zero-condition" is equivalent to the so-called "sum-rule" and close to 
the Strang-Fix condition of the same order. It will be convenient for us to 
define Z~ x as the set of all trigonometric polynomials. 

3. Auxiliary results 

Proposition 1 A trigonometric polynomial t belongs to Z° if and only if 

t„(0) = ^, u = 0,...,m-l, (3) 

m 

where r , . . . , r m _i are the polyphase functions oft. 
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Proof. Let s £ D(M*), using (2) we have 

m— 1 m— 1 

t(M*~ 1 s) = J2 e2m{r "' S) Ms) = ^e 2 ^ r - s V,(0). 
So, the relation t G Z° can be rewritten as 

m— 1 

]T e 2m ^r u (0) = t(0)5 0l , 1 = 0, . . . , m - 1. 

Consider these equalities as a linear system with unknowns To(0), . . . , r m _i(0). 
Due to Proposition A, the system has a unique solution given by (3). 

Lemma 2 Let t be a trigonometric polynomial, t(0) = 0, then 

d 

t(x) = J2tk(x)(l-e 2m ^), 
k=i 

where t k , k = 1, . . . , d, are trigonometric polynomials. 

Proof. There exists N e Z d such that t(x)e 27ri{N ' x) = p(z), where z = 
exp(2nix) and p is an algebraic polynomial. It is clear that P(l, . . . , 1) = 0. 
By Taylor formula it follows that 

d 

P( Z ) = ^2Pk(z)(l - Z k ). 
k=l 

where p k , k — 1, . . . , d, are algebraic polynomials. It remains to set t k (x) : = 
p k (z)e- 2 ^ N ' x \<) 

Lemma 3 Let t,t G Z°, t(0) = i(0), then 

d 

t(x) - t(x) = JTt k (x) (1 - e 2m{M * x ^) , (4) 

k=l 

where t k , k — 1, . . . , d, are trigonometric polynomials. 
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Proof. Using (2) for t, t, we have 

m— 1 

t(x) - t(x) = e 2 ^ x \r v {M*x) - f v {M*x)). (5) 

Due to Proposition 1, 7^(0) — t„(0) = 0, f = l,...,m — 1. Hence, by 
Lemma 2, 

r„(y) - Uv) = E ( x " e2m( "' efc) ) ' ( 6 ) 
i=i 

where T„k, k = l,...,d, are trigonometric polynomials. It remains to set 
y = M*x and combine (6) with (5). 

Corollary 4 Let t E Z°, r EZ d ,t := e 27ri ^H, then (4) holds. 

Lemma 5 Let t E Z~ x , r E Z d , t := e 2nt<yT ''h, then the k-th polyphase 
function oft is given by fk(x) := e 2m ^ nv ' x ^T v (x), where v is the number of the 
unique digit s u E D(M) such that s u + r = Sk + Mn v , n v E Z d . 

Proof. By (2), 

m—l m—1 

t(x) = e 2wi{r ' x) e 2 ™ {Sv ' x) T v (M*x) = ^ e 2wi{s " +r ' x) T u (M*x). 

It is clear that s v + r — s tv + Mn u , n u E Z d , {/ , • • • , l m -i} = {0, . . . , m — 1}. 
So, 

m—l 

t( x ) = Y,e 2m{siv ' x) e 2ni{n '" M * x) Tv{M*x), 

and the trigonometric polynomial fi v {x) := e 2mi ^ lv ^ r v {x) is the ^-th polyphase 
function of t. 

Lemma 6 Let n E Z d , 1 - e 2ni{M ' ln ^ ] = for each j = l,...,d. Then 
n = (mod M). 

Proof. Since e 2 ™ = 1 if and only if n E Z, we have (M~~ 1 n)j = 
(M~V ej) E Z for each j = l,...,d. Hence, M~ l n = lEZ d , i.e. n — Ml. 
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4. Decomposition of masks 

In this section we prove the existence of a decomposition of a trigonometric 
polynomial t G Z n of the form 

d 

(1 - e 2 ^ x ' ek) )t{x) = ^2t jk (x)(l - e 2 ^ M * x ^), 
3=1 

with tjk G Z" 1 ' 1 ,j = 1, . . . , d, for all k = 1, . . . , d, and present algorithms 
for its computation. We start with the simpler case n = 0, 1, and then do 
the general case n > 1, which is much more complicated. 

4.1 The case te Z n , n = 0, 1 
Proposition 7 Lei t G £/ien 

(1 _ e 2ni( X ,e k )^ = J2t jk ( X )(l - e 2 ™( M *^)), (7) 
3=1 

where tj k , j = l,...,d, are trigonometric polynomials. Conversely, if (7) 
holds for a trigonometric polynomial t, thent G Z° andtj k (0) = (M~ r )j k t(0). 

Proof. If t G Z°, then (7) follows immediately form Corollary 4. Now let 
t be an arbitrary trigonometric polynomial such that (7) holds, s G D(M*), 
s 7^ 0. It follows from (7) that 

d 

t{M*~\s) (l - e^m*- 1 ^))^ = ^t jfc (M*- 1 s)(l - e 2 ™^)) = 0. 

3=1 

Taking into account Lemma 6, we obtain t(M*~ 1 s) = 0, which means that 

tez°. 

Rewrite identity (7) as follows 

d 

(1 - e 2m{x > M ~ lek) )t{M*~ l x) = ^t^(M*- 1 x)(l - e 2 ™(^')), 

3=1 
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Differentiating by x h we have 

-2m (M-^fe, e i )e 2m ^ M_lefc ))t(M*~ 1 x) + (1 - e 2m{x ' M ~ lek) )-^-t{M*- l x) = 

OXi 

d Q 

^ -^t jk {M*- l x){l - e 2 ™^) - 2int lk (M*- 1 x) e 2m{x > ei) . 
j'=i ' 

Substituting x = 0, we obtain 

^^(M-^e^O) = -27Tit, fc (0). 

It remains to note that (M _1 e fc ,e ; ) = (M~ v ) lk .() 

Analyzing the proofs of Lemmas 3 and 5 it is not difficult to describe an 
algorithm for finding the polynomials tj k , j,k = 1, . . . , d, in the decomposi- 
tion (7). We assume that the polyphase functions r v of t are given, and the 
algorithm derives the polyphase functions Tj ku of each polynomial tj k . 

ALGORITHM 1 

Input: {?„(/), leZ d , v = 0, . . . , m - 1}. 

Output: {jjkvip)i I ^ u — 0, . . . ,m — 1, j, k — 1, . . . , d}. 

Step 0. For each v — 0, . . . , m — 1 find the sets 

n v := {leZ d : + 0}. 

For each k — 1 , . . . , d do 
For each v — 0, . . . , m — 1 do 

Step 1. Compute M' 1 ^ — s u + s n ) =: Z„ for all n — 0, . . . , m — 1 
and denote by n* the unique n such that Z„ G 7h d . 

Step 2. Find the set fl ku :— {I — r + l n *,r G fi n *}, and for each 
/ G put f„(Z) := ?„»(/ - /„*), 

Pkv{z):= ^ (f„(Z) -f v {l))z l . 
ien v un kl , 

Step 3. Put 

Tlkv(x) = — (Pku(z) - Pfc„(l, Z2, ■ ■ ■ , , 

1 — Zi 

T2kv(x) = ~, (Pfci/(1, *2, Z3, . . . , Z d ) - p fc „(l, 1, Z 3 , • • • , Z d )) , 





1 


1 


- Zi 




1 


1 


- z 2 




1 


1 


- 2d 



Tdkv(x) = z (Pfci/(1, 1, • • • , 1, Zd) - Pfci/(1, 1, • • • , 1, !)) -0 



Theorem 8 If t G Z 1 , then in any decomposition (7) the trigonometric 
polynomials tj k , j, k — 1, . . . , d, belong to Z°. 

Proof. Let I = 1, . . . , d, s G D(M*), s ^ 0. Since t G Z {1 \ 
d 



V_ ^ _ e 2ni(M*-i X ,e k ))^ t ( M *-l^ 



0, k=l,...,d. (8) 



On the other hand, it follows from (7) that for each k — 1, . . . , d we have 

_d_ 

dxi 



dxi 

d 

*-l 



E^( M *" ls )^( 1 - e27r ^) 



= —2mtik{M*~ L s). 



So, we proved that ti k (M*~ 1 s) = for all s G D(M*), s ^ 0, which means 
4.2 The general case t G Z n , n > 1 

Proposition 7 and Theorem 8 state that for n — 0, 1 the condition f G 2™ 
implies tjfc G -Z" -1 . This fact can not be extended to the case n > 1. 
If t G Z 2 , there exist decompositions (7) whose elements are not in Z 1 . 
Nevertheless, it will be shown that any decomposition of t G Z n can be fixed 
to provide t jk G Z^ 1 . 

Theorem 9 Let n G Z + . A trigonometric polynomial t belongs to Z n if and 
only if there exists a decomposition (7) with tj k G Z n ~ l , j,k = 1, . . . , d. 

Proof. Set 

a(x) = t(M*~ l x), a jk (x) = t jk (M*~ 1 x), 
b k (x) = l- e 2 ^ M *~ lx ^ , c k (x) = 1 - e 2m{x ' e ^ . 

In these notation, (7) can be rewritten as 

d 

b k (x)a(x) = ^a jk (x)cj(x). (9) 

3=1 



9 



Note the following trivial properties of cf 



Cj (l) = 0Vl£ Z d , (10) 
Cj = Vc> 7^ rej, 



D 5 Cj = V<5 ^ re,, r G N, 5 + 0, (11) 



p-(l) + V/ G Z d . (12) 

ax,,- 

We will prove the theorem by induction on n. We have a base for n = due 
to Proposition 7. Let us prove the inductive step: n — > n + 1. 

Assume that each polynomial belongs to Z n , i.e. D^ajk(sl) = 0, 
Z = 1, . . . ,m - 1, for all (3 G Z d + , [0\ < n. Let a G Z d , [a] = ra + 1, 
s G D(M*), s^O. It follows from (9) and Leibniz formula that 



E " 

0</3<a 



) D a -^ fe ( S )^a( S ) = D a a jfe (s) Cj ( S ). 



By the inductive hypotheses, D^a{s) = whenever < a. Hence, taking 
into account (10), we have 

b k (s)D a a(s) = c 3 {s)D a a jk {s) = 0. 

Due to Lemma 6, b k {s) ^ for at least one k — 1, . . . ,d. It follows that 
D a a(s) = 0. 

Now we assume that t G Z n+1 , i.e. D a a(s* v ) = 0, v = l,...,m — 1, 
for all a G Zl, [a] < n + 1. Due to the inductive hypotheses, there exist 
polynomials tj k G Z n ~ l satisfying (7), i.e. D 5 aj k {s* p ) = 0, v — 1, . . . ,m — 1, 
for all 5 G Zl, [5] < n — 1. We will construct new trigonometric polynomials 
satisfying (7) and belonging to Z n . Note that for n — any trigonometric 
polynomials satisfying (7) belong to Z°, due to Theorem 8. 

Let us prove the following statement. If D 5 aj k {s* v ) — 0, v — 1, . . . , m — 1, 
for all 5 G Z+, [5] < n and for all j = 1, — 1 (/ = 1, . . . , d), then there 

d d 

exist polynomials tj k , j,k = 1, . . . , d, such that Yl ®jk c j — Yl a jkCj, where 

3=1 3=1 

dj k (x) = t jk (M*~ 1 x), and D s a jk (s* u ) = 0, v — 1, . . . , m - 1, k— 1, . . . , d, for 
all 5 G , [5] < n, and for all j = 1, ... , /. 
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Let s G D{M*),s ^ 0, (3 G Z d + , [/?] = n, £ # = 0. Set a = (3 + e t and 

i=i+i 

note that [a] — n + 1. It follows from (9) that 



X 



D a (b k (x)a(x) 



0. 



(13) 



On the other hand, due to (10), (11) and the assumption of the statement, 
we have 



^D^a jk (s)^(s) = D% k (s)^(s) 



dxi 



Combining this with (12) and (13), we get D"ai k {s) = 0. is proved already 
for I — d (in this case tj k = tj k , j,k — 1, . . . , d). 

d 

Next let j = I + 1, . . .,d, (3 e 1%, [(3} = n, £ # > 0. Define the 

i=i+l 

functions 

j m— 1 

5/Jj(aO = 9n-i,p- Bj (x) h u {x)D l5 ai k (sl), 

u=l 



where g^s is a trigonometric polynomial such that Vg^siO) = for all 
7 G , [ 7 ] < N, 7 ^ 5 and D s g m (0) = 1; 



1 m— 1 

h< X ) = - V e 2^(-<,M-is M ) 



Since, by Proposition A, h„(s*) = S^, due to (10), (11) and Leibniz formula, 
it is not difficult to see that for each v — 1, ... ,m — 1 we have 



D^ Cj {x)q Pj {x)) 
D 5 ( Cj (x)q Pj (x)) 
D s (ci(x)qpj(x)) 



e Z d + , [5] <n, 8^(3; 



(14) 
(15) 



X s. . 



V5eZ d + , [5] < n- l,Vi= (16) 



Set 



a ifc (x) := a ifc (x) - ^ c j (x)g /3j (x), 



j=;+l [/3]=n, 0j>O 

l+1 = ...=l3 j _ 1 =O 
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a jk (x) := a jk (x) + £ Q(ar)g^j(x), j = Z + 1, . . . , d. 

[0]=n, /3^>0 
0J+l=-=0j-l=O 

d d 

Because of construction, E djk c j = E a jfc c j> an d, taking into account (14), 
(15), for each = 1, . . . , m — 1 we have D^ai k {s* v ) = whenever [/?] = n, 
E A > 0; Df>a lk (s* v ) = D^a lk (s* v ) = whenever [/?] = n, £ # = 0. At 

i=H-l i=J+l 

last, due to (16), D s a jk (s* u ) = 0, j = /,..., d for all 5 G Z*, [<J] < w - 1. 
To complete the proof of the statement it remains to put tj k = dj k (M*x) for 
j = I,. . . ,d and t jk = t jk for j = 1, ...,/- 1. 

So, we described one step for improvement of decomposition (7). Starting 
with 1 = 1, after (d — 1) steps we will obtain required polynomials. 

Analyzing the proof of Theorem 9 it is not difficult to describe an al- 
gorithm for finding polynomials tj k G Z n ~ l , j,k = 1, . . . , d, in decomposi- 
tion (7). To realize this algorithm we will need the functions g^s- Explicit 
recursive formulas for these functions are presented in [10]. 

ALGORITHM 2 

Input: t G Z N , N > 1. 
Output: tj k G Z^ 1 , j, k — 1, . . . , d. 

Step 1. Using Algorithm 1, find tj k , j, k — 1, . . . , d, satisfying (7). 
Step 2. For n = 1, . . ., N- 1 do 

Set tfk := tjk, j,k = l,...,d; 
For each I — 1, ... ,d — 1 do 
For each k — 1 , . . . , d do 

a lk (x) = tt 1 \M*- 1 x); 

j=l+l [0]=n, 0j>O 

For each j — 1, . . . , I do 



For each j — I + 1, . . . , d do 
== + (1 - e^"**^) £ ^(M^).O 



[/3]=n, /3^>0 
^+1=-=^_1 = 
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In the case d = 2, Step 2 of Algorithm 2 does not look so frightening, it 
is reduced to the following. 
For n — 1, . . . , N — 1 do 
Set 1 = 1, t- k '.= tjk, j,k = 1,2; 
For each k — 1, 2 do 



t lk (x) := ^(x)-(l-e 2 ^-)) J2 ^( M **); 

[/3]=n, /3 2 >0 



=n, /3 2 >0 



In particular, for d = 2, iV = 2, Step 2 may be realized as follows. 

t/fc := tjk, j, k — 1, 2; 
For each = 1, 2 do 



m—l 



t lk (x) := 4l\x) + ± (1 - e 2 ^*) £ ^(M*x) A t (0) (M -i u) 

1 m — 1 rj 

t 2fc (x) := - ^ (1 - e^'^) £ K(M*x)^(M*-\) 

For each n G N we introduce the set 

r n := {fc G I" : jfc, g {1, . . . , d}, l = l,...,n}. 

Theorem 10 Let n, n G N, n < n , t G Z n °~ l ,then there exist trigonomet- 
ric polynomials tjk G Z n °~ n ~ x , k,j G T n , such that 

n d d n 

H(l-e^^)t(x) = ■ ■ ■ E II C 1 - e 2 ^ M ^) , (17) 

J=l J'l=l J'n=l i=l 

n 

t iJfc (o) = n( M-i u*(°)- 
i=i 
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Proof. We will proof by induction on n. 

Base: n — 1. Let 6 T 1 . Due to Proposition 7 and Theorem 9, there 
exist trigonometric polynomials tj k G Z n °~ 2 , j G T 1 , such that (7) holds. So, 
we have (17) for n — 1, (18) also follows from Proposition 7. 

Inductive step: n — 1 — > n. Let 1 < n < n , A; G T n , fc' := (fci, . . . , &„_i). 
By the inductive hypotheses there exist tf k > G Z n °~ n , j' G r™ _1 , such that 

n— 1 (i d n— 1 

J-J ^ _ ^(s.e*,)) ^ = t JJ (1 _ e 2-(M^,e 3i )^ ^ (lg) 

1=1 jl = l J'n-l=l J=l 

n-1 

wo)=n( M-i u*(°)- ( 2 °) 

Since n — n > and the theorem is proved already for n — 1, for each fy fe /, 
there exist trigonometric polynomials t jnkn G .g^o-n-i^ j e pi ; suc h that 

(l-e^O)^) = ^^(x)(l-e 2 ^ M *^)), (21) 

jn=l 

^fen(O) = (M- 1 ) jnfc „V fc ,(0). (22) 

Combining (19), (20) with (21), (22) we comlete the proof. 

To impart a more compact form to (17), (18), we introduce the following 
notations. Set 

A k (x) = A k = (1 - e **(e*.*)) , A(x) = A = (Ax, . . . , A d ) T , 
5 k (x) = 5 k = (1 - e 27r ^ M **)) , =* = (*!,..., fcf. 

Now Theorem 10 can be rewritten as 
Theorem 10' Let n, n G N 7 n < n , t G Z n °~ l ,then there exists a d n x d n 
matrix T whose entries are trigonometric polynomials T k j G Z no ~ n ~ 1 , for 
k,j G T n , such that 

(A(x)) [n] t(x) = T(x)(5(x)) [n] , T(0) = t(0)(M*- 1 ) [n] . 
5. Construction of masks in Z n 

A simple description of the classes Z n is well known in the one-dimensional 
dyadic case. A general form is given by the formula t(x) = (1 + e 2mx ) n T(x), 
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where T is an arbitrary trigonometric polynomial. In the multidimensional 
case Z n can not be described in a similar way. We will give a characterization 
of the class Z n for arbitrary dilation matrix which allows to construct its 
elements in practice. 

Theorem 11 A trigonometric polynomial t belongs to Z n if and only if the 
derivatives of its polyphase function r k , k — 0, . . . , m — 1, up to order n are 
given by 



£><V fc (0) = - J2 xJ a g )(-2mr k r-^ a E Z+, [a] < n, 

^ n<- R<r n < \P / 



(23) 



0<f3<a 

where X a = D a t(M*~ 1 x)\ x=0 . 



Remark. For polynomials t whose polyphase functions To, ... , r m _i form 
a unimodular row (i.e. there exists a dual row of trigonometric polynomials 
To, ... , T m _! such that YlT^o T kTk = 1), the statement of Theorem 11 follows 
from combining the results of [2] and [10]. It was proved in these papers that 
both conditions are equivalent to vanishing moments of the corresponding 
wavelet system. 

Proof. Assume that (23) holds with some complex numbers A a . Let 
s G D(M*). By Leibniz formula, 



D a (e 2m{rk ' x) r k (x)) = V ( a \ (e 27Ti{rk > x) ) D a ^r k (0) = 

X=S 0</3<cA/V 

(f) e 2m{r *< s \2mr k yD a -?T k (0) = 



0</3<a 



E (fji^f E A 7 ( a :") ( -2„,,r^ = 

e e A T ( a -^(;)(-2„,,)-n(-ir^ 

<0<aO<y<a-0 V ' / \ H / ?= l 



_ 2wi(r k ,s) 

m 

O<0<a x ' ' 0<7<a-/3 

J_ e 27ri(r fe ,s) 

m , , 

()• I- a ()• a i j 

0<7<a V ' 7 0</3<ct-7 \ ^ / j=l 
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Since 



E 

0</3<a-7 



a — 7 

p 



n(-i)- ft 



n E 

j=l 0<Pj<aj—jj 

f[(i - ir 



A' 



0, a 7^7, 

1, a = 7, 



we have 



D a (e^^V^x)) 



A 



« e 2«(r t ,.) j fc = o,..., m -l. 



It follows from (2) and Proposition A that 

m—l 

L> Q (t(M* _1 :r)| = D a (e 27ri{rk ' x) r k (x)) 



k=0 



m 



m—l 

E' 

fc=0 



,27ri(r fc ,s) 



A Q , if s = 0, 
0, if s^O. 



Now let us check that (23) follows from the relation t G Z n . We will prove 
by induction on n. The base for n = was established in Proposition 1. To 
prove the inductive step n — > n+ 1, we assume that £ G 2 n+1 and (23) holds. 
Let aeZ d , [a] = n + 1, s G D(M*). By (2) and Leibniz formula, 



E E ("jD^(e 2 ^^)| D"t*(0) = 

m—l / 

E ( + E 

0</3<a 



D Q (£(M* x) 



k=0 



D a - l3 (2nir k ) a - f3 D /3 r k (0) e 2m{r ^ s) . (24) 



Because of the inductive hypotheses, we have 



m E E X ~ 

0</3<o 0<7</3 



E 

0</3<o 

cA / (3 



(2mr k ) a ^D^r k (0) = 



PJ V7 



(2mr k ) a ^(-2inr k ) 
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d 



;E^rs6)(5n<- 

0<7<q 7</3<a v ^ 7 v ' 7 j=l 

; e (;) w e (v)n<-^ 
k e (:)A 7 (2^r-(ri(i-ir— -n(-if) 

0<7<a V ' 7 \j"=l j"=l / 



i e (;) 

0<7<o v ' 7 

Combining this with (24), we obtain 

J2e 2m ^ s UD a T k (0)-^ J2 I ) A 7 (-2^r fc ) Q - 7 ) = D a {t{M*~ l x)) 

fc=0 \ m 0<7<o^^ 7 / 

Set A a = Z} a (i(M* -1 a;))| x=0 . Due to Proposition A, the linear system 

m— 1 

£ e 2«(r*,.r) yfc = A a 5 0i , / = 0,...,m-l, 



fc=0 



has a unique solution y& = fc = 0, . . . , m — 1. It follows that 

^ ^ j A 7 (-27rzr fe r~' = >V 

0<7<a 



m ^1 \ ry I 



m 



which was to be proved. 

So, if we want a polynomial t to belong to Z n , its polyphase functions 
should have derivatives at the origin given by (23). This can be easily realized 
for an arbitrary set of parameters \p, [f3] < n. General forms for all such 
polynomials t are presented in [10]. 



6. Subdivision schemes 

In this section we apply the decompositions of Section 4 to the analysis of 
convergence and smoothness of multivariate subdivision schemes associated 
with general dilation matrices. We first present basic definitions, and prove 
an important observation about the matrices in a decomposition of t G Z 1 . 
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6.1 Preliminaries 

Let A a , a G Z d , be N x N' matrices such that A a is non-zero only for a finite 
number of a, then T{x) := Yla& d A a e~ 2nt( - a ^ is & N x N' matrix whose 
entries are trigonometric polynomials (N x N' trigonometric matrix). The 
subdivision operator St = <St,m associated with T and with a dx d dilation 
matrix M is defined on £^ = £^ (Z d ) by 

(Srf)a = ^2 Aa-Mpfp, f £ ^o- 

It is clear that St is a linear bounded operator taking £^ into l^. If N = N', 
the operators S%, n = 1, 2, . . ., are well defined, and the sequence {S T }%L 1 is 
& N x N matrix subdivision scheme (scalar subdivision scheme if N = 1). 

Let t be a trigonometric polynomial, T be a d x d trigonometric matrix, 
set L(x) = (Li(x),...,L d (x)) T , where L k (x) = (1 - e 27ri{x > ek) )t(x), R(x) = 

(Ri(x), R d (x)) T , where R k {x) = £ T kj (x)(l - e 2 ^ M ' x ^). To each / G 

j'=i 

£oo assign the vector-valued sequence V/ G ^ defined by (V/) a = (/ a — 
/a-ei, ■ ■ • , /a - fa- ed ) T ■ It is clear that 5 L / = V ,%(/), S^/ = 5 T (V/) for 
any / G 4x>- Hence equality (7), with = Tj^, may be rewritten in the form 

v5 t (/) = 5 T (v/), V/e4o. (25) 

Similarly, to each / G ^ we assign V / G defined by 

(V/) a =((v(/,e 1 ))^...,V(/,e JV ))a T . 

If all trigonometric polynomials tj k in the right hand side of (7) are in 
Z°, we can decompose them (see Theorem 10). This second step of the 
decomposition may be rewritten as 

VS T (f) = S Q (Vf), V/GC (26) 

where Q is a d 2 x d 2 trigonometric matrix. 

Proposition 12 Let t be a trigonometric polynomial, t G Z 1 , t(0) = m, and 

let T be a dxd trigonometric matrix satisfying (25), T(x) := X]«ez d A a e~ 2m ( a ' x \ 
Then 

A a - Mf 3 = M*" 1 \faeZ d . 
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Proof. First of all, we note that 2^/3 e z d A*-M/3 = X^e^ A*'-m/3 when- 
ever a = a' (mod M). So, it suffices to check that 

X u :=Y, A s,-M P = M*-\ i/ = 0,...,m-l. 
Substituting x = M*~ 1 s* k into the equality 

m—1 

T(x) = Y^e- 2m ^ A Sv - M pe- 2 ^ M ^ x \ 

we have T(M*" 1 s*) = YZ=o e 2 ^^' 1 ^ X v It follows from Proposition 7 
and Theorem 8 that 

m—1 

e -2^,M"-'si) Xv = m s k0 M*-\ k = 0, . . . , m - 1. 

This linear system has a unique solution X„ = M* _1 , z/ = 0, . . .,m — 1, 
because of Proposition A.O 



6.2 Convergence 

A scalar subdivision scheme St = St t M, associated with a dilation matrix M, 
is called uniformly convergent, if for any / G £oo, there exists a continuous 
function S^ 00 / such that 

lim ||^/-Sr/(M- fc .)||oc = (27 ) 

fc— >oo 

and if for at least one / G £oo, the limit function S%°f is not identically zero. 

The issue of the convergence of a multivariate scalar subdivision scheme 
associated with a dilation matrix M is studied in [4]. The convergence 
result there is limited to dilation matrices having a self-similar tile. A d x d 
dilation matrix M is said to have a self-similar tile if there exist a set of digits 
D(M) and a bounded set E C M d whose integer translates form a disjoint 
decomposition of lR d such that 

ME = U(E + s), (28) 

seD(M) 
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Theorem 13 [4] Let S t = S^m be a scalar subdivision scheme and let there 
exists a self-similar tile related to M . Then S t is uniformly convergent if and 
only if: a) t(0) = m; b) there exists a d x d trigonometric matrix T such 
that (25) holds, c) lim sup \\S%S! f\\ id =0. 

fc^oo feloo °° 
l|V/||oo=l 

Here we present sufficient conditions for convergence of a scalar multi- 
variate subdivision scheme associated with any dilation matrix M, and then 
give an example of a subdivision scheme satisfying these conditions, where 
the verefication of the conditions is done with Algorithm 1. 

Theorem 13' Let t G Z° and M be a d x d dilation matrix. Then there 
exists a trigonometric matrix T of order d x d satisfying (25). Moreover, if 
t(0) = m and 

lim Halloo = 0, (29) 

then the subdivision scheme S t is uniformly convergent. 
First we prove a simple lemma. 

Lemma 14 Let n G Z +; M be a d x d dilation matrix, and let ip be a 
compactly supported function satisfying the refinement equation 

<p(x) = a a (f{Mx - a), Vi 6 K d . (30) 

a<=Z d 

Then 

Y,fMx-P) = Y,( S *nf)MMx-a), VxGR rf (31) 
for any f G t^, where 

i(x) = a a e~ 2m{x ' a) 

Proof. Using (30), we have 

(Si In f) a <p(Mx - a) = Y Y a «-Mf]fp ¥(Mx - a) = 

/3ez d aez d /3ez d 
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Proof of Theorem 13'. The existence of T satisfying (25) follows from 
Proposition 7, as indicated in Subsection 6.1. 

Let <f be a continuous compactly supported function satisfying (30) and 
the interpolatory conditions 

<p{a) = 6 a0 , WaeZ d , (32) 

The existence of such a function for an arbitrary dilation matrix M is proved 
in [3, Proposition 4.1]. It is proved in [3] that the corresponding mask t is in 
Z° and that t(0) = m. To prove the second part of the claim, we will show 
that the sequence of functions 

F k (x) = ^(M fe x-a)(S*/)a, fc = l,2,..., 

aeZ d 

is a Cauchy sequence in C(lR d ) for any / G £oo- 
By Lemma 14 with n — 1, we get for / G £oo 

F k (x) = <P(M k x - a)(S*f) a = ^M k+1 x - a)(S~ t S*f) a , 

at=1 d a<=1 d 

and therefore 

F k+1 (x) - F k {x) = ^{M k+1 x - a)((S t - 5 f )5*/) a . (33) 
Since t and t are in Z° and t(0) = t(0), we get from Lemma 3 that 

d 

t(x) - t(x) = Qk(x) (1 - e 2 ^ M ^) . 
k=i 

This leads, by similar derivations to those leading to (25), to the existence 
of a vector of subdivision schemes (S qi , . . . , S qd ), such that for any / G £oo, 

(S t -S~ t )f = S t _~J=(S qi ,...,S qd )Vf. 

Hence, (33) can be rewritten as 

F k+1 (x) - F k (x) = J2 i(S qi , S qd )vS*f)MM k+1 x - a), 
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Using relation (25), we obtain 



\F k+1 (x) - F k (x) 



((S qi , • • • , S qd )S k T Vf)MM k+1 x - a) 



< 



qi^v/lloo^cn^iuiv/i, 



where C depends on ip, q ± , . . . , q d . 

Condition (29) is equivalent to the existence of a positive integer L such 
that H-Srlloo = fJt < 1. Thus 

\F k+1 (x) - F k (x)\ <C\7f\^. 

This yields that for all k, n e N 



n-l 



\F k+n (x) - F k (x)\ < C\Vf\J2^ < T^IV/I^, (34) 



which implies that {F k } is a Cauchy sequence. Denote the limit function by 
F which is, evidently, a continuous function. Passing to the limit in (34) as 
n — > oo, we have 

(35) 



lim sup \F(x) - F k (x)\ = 0, 



Due to (35), 



lim sup 

fc^oo 



F(x)- J2(Stf)MM k x-a) 



Substituting x = M k f3, (5 £ Z d , we have 



lim sup 

k^oo 



0, 



The claim of the theorem follows now from (32). 

2 
2 -1 



Example. Let M 



,s 2 







,s 3 
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Define the polyphase functions of a polynomial t by 



t (x) = — (A + Az 1 +Az 2 + 4z 1 z 2 ), 



Tl(x) 
T2(x) 



16 
1 

16 
1 



(5 + 4*i + z{ x + 2z 2 + 3z 2 x + Zl z 2 ), 
(4 + z x + 2zf 1 + 5^2 + z 2 _1 + 3-2i z 2 ), 



t 3 (x) = -^(5 + zx + Az 1 1 + z 2 + 3z 2 x + z 1 z 2 + z 1 1 z 2 ). 



It is not 
we find 


difficult to see that t £ Z°, and that t(0) = m 


Using Algorithm 1, 


Tuo(x) 


= ^(-1 + 2^2 + z x z 2 + 2z 2 2 + Zl 4), 


T 2 lo(x) 


= ^(5 + 3z 2 ), 


t u1 (x) 


= Te( z i 1 + 3z 2 ), 


T21l(x) 


= 1^(5 + 3z 2 1 ), 


Tu 2 (x) 


= ^(-1 + 2z^ + 2 2 + 2 2 2 - z^z 2 ), 


r 2 i 2 (x) 


= ±(5 + 3z 2 + z 2 1 ), 




= ±(2z^ + 2z 2 + z^z 2 ), 


r 2 i 3 (x) 


= ±(5 + 2z 2 l ), 


Tl2o(x) 


= + zi + 4z 2 + z 2 x + 3ziz 2 ), 


r 22 o(x) 


— 16 Z 2 ; 


Tl2l(x) 


= jq(2 + z 1 + z{ x + z 2 + 3z 2 x + ziz 2 ), 


T 22 l(x) 


= 0, 


n 22 (x) 


= ±(3 + Z2 + 2z; 1 ), 


T 222 (x) 


— 16^2 > 




= ±(3 + 3z^ + z^z 2 ), 


t 223 (x) 


= 0. 



By Theorem 13', to prove that the subdivision scheme S t is convergent, it 
remains to check that the norm of the matrix subdivision operator St is 
strictly less than 1. This requirement is fulfilled because 

15 

H-StIU < maxdlrifc^H^ + ||r 2 ^||^) < — . 

(u,k) 10 

(Here we identify a trigonometric polynomial with its sequence of Fourier 
coefficients) . 
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6.3 Smoothness 

Now we discuss how to study the smoothness of the limit function of a uni- 
formly convergent scalar subdivision scheme. 

Let t — a a e~ 2m ( a ^ be a trigonometric polynomial, and let t„, v = 

aez d 

0, . . . , m — 1, be its polyphase functions. It is well known (see, e.g., [4]) that 

r ^(°) = a ^~ M P =1 ' ^ = 0, . . . , m - 1, 

/3ez d 

whenever (27) is fulfilled for at least one / G £oo f° r which S^ 00 / ^ 0. So, if 
St uniformly converges then, due to Propositions 1 and 7, there exists a dx d 
trigonometric matrix T such that (25) holds. 

Following [9] , we introduce the notions of a normalized subdivision scheme 
and its subconvergence. 

For X a NxN matrix, and T a NxN trigonometric matrix, the sequence 
{X* k S T } ( £ =1 is called the normalized (by X) subdivision scheme St- 

We say that a subdivision scheme St = St,m normalized by X is uni- 
formly convergent on a subspace L of L ^ {0}, if for any F <E L, there 
exists a continuous vector-valued function S^, X F (S T °, x F(x) G M^) such 
that 

Urn \\X* k S k F - S T ° /x F{M- k -)\\ OQ = 0. (36) 

and if for at least one F G L the limit S T /xF is not identically zero. 

We say that a subdivision scheme SV normalized by X is uniformly sub- 
convergent on a subspace L, L C L 7^ {0}, if for some infinite set C N 
and for any F G L there exists a continuous vector-valued function S^-^F 
(S%> /x F(x) G R d ) such that 

lim ||M* fc ^F - S^ /x F(M- k -)\\ 00 = 0. (37) 

and if for at least one F G L, the limit St/xF is not identically zero. 

Next, following [8], we consider the class of isotropic dilation matri- 
ces. A matrix M is called isotropic if there exists a constant C such that 
||M fc || 0O ||M- fc || 0O < C for all k G N. 

Here we give a detailed proof of the sufficiency of a necessary and sufficirnt 
condition for the limits of a scalar subdivision scheme to be in C 1 . This 
condition is stated in [9] with a sketch of a proof, which is not applicable 
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for general isotropic dilation matrices. Indeed, by Proposition 12 here, St 
cannot be of full rank, as required by the sketch of the proof, whenever M 
is not a multiple of the identity matrix. 

Theorem 15 Let M be isotropic, let a scalar subdivision scheme St = S^m 
be uniformly convergent, and let St be a d x d matrix subdivision scheme 
satisfying (25). If the subdivision scheme St normalized by M uniformly 
subconverges on V^, then the function f is in C 1 (M d ) for all f G £oo- 

First we prove a lemma. 

Lemma 16 Let 9JT be an infinite subset of N, g k G P^, k G 9JT, g be a 
continuous vector-valued function (g(x) G M. n ), and let fy k , k G 9JT, be n x n 
matrix-valued functions which are uniformly bounded, uniformly compactly 
supported and such that 

J2^ k (--a) = I n , VkeWl. (38) 
lim \\g k -g(M~ k -)\\ oo = 0, (39) 

k— >oo 

then the sequence of vector-valued functions J2 a <=z d ^k(M k -—a)(gk) a , k G 971, 
uniformly converges to g on any compact set K C M. d . 

Proof. Let K G M. d be a compact set, x G K, k G 9Jt. Because of (38), 
^k(M k x - a)(g k ) a - g(x) = £ V k (M k x - a) {{g k ) a - g(x)) = 

V k (M k x - a) {{g k ) a - g(M~ k a)) + 

a<=n(M k x) 

^k(M k x - a) (g(M~ k a) - g(x)) , (40) 

where Q(t) = U Sl k (t), with Sl k (t) = {a G Z d : tf fc (t - a) ^ O n }. 

It is clear, that ft < C, where C is a constant depending only on the 
joint support of the functions ^ k . It follows from (39) that 



lim sup 

K- — * OO 



^t{M k x - a) ((g k ) a - g{M~ k a)) 



0. (41) 



25 



If q k (M k x -a) ^ O n , then \M~ k a - x\ < \\M- k \\R, where R is the radius 
of a ball containing the supports of all Due to the uniform continuity of 
g on a compact set, this yields 




To complete the proof it remains to combine (41) and (42) with (40). 
Proof of Theorem 15. For a € N d , the function 

d 

b A x ) = Y[(X[o,i] * •••*X[o,i])(xj) 

is the tensor product B-spline of order a, when the number of convolutions in 
the above product are Oj + 1. In the following we take <jj > 1, j — 1, . . . , d. 
Thus, b a has the following properties: 

b*(x - a) = 1, (43) 

a&L d 



grad b a = (b a - ei (• - e x ) - 6 CT _ ei , . . . , 6 <y _ ed (- - e d ) - & CT -e d ) • 



(44) 



Given f E L, using (43) and Lemma 16 with 9JI = N, n — 1, ^ = S k f, 
g = S%°f, we obtain that the sequence F k := J2 a ez n ba(M k ■ -a)(S k f) a 
uniformly converges to S^°f on any compact set K C W 1 . These functions 
are m C l {M. d ) by the choice of a. To prove the theorem it remains to check 
that the sequence of vector- valued functions grad F k uniformly converges to 
a continuous vector-valued function on any compact set K C W 1 . 
Set 



( k-ei 



V 



o 



b a -e d 



Using (25), we have 

grad J2 h ^ M " x - a )( S tf)a = M* k J2 g™db a (M k x - a)(S k f) a = 

M* k B a {M k x - a)(vS k f) a = M* k B a (M k x - a)M*- k (M* k S£v f) a . 
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Due to the normalized uniform sub convergence of St on VL, there exists 
an infinite set Dt C N such that for any F G VL there exists a continu- 
ous vector-valued function S^^F for which (37) holds. Let / G L. Since 
Eae^ M* k B a (x - a)M*~ k = I d , by Lemma 16 with JOT = 0t, n = d, 
gk = M* k S T yf, g = 5~ /M V/ and ^ k = M* k B a M*~ k (the functions # fe 
are uniformly bounded because M is isotropic), we can state that the se- 
quence gradF fc , k G 91, uniformly converges to S^ M \7f on any compact set 
K cR d . 

A direct conclusion from Theorem 15 in case equation (25) holds for 
a trigonometric polynomial t, is that the limit function of the subdivision 
scheme S t is in C 1 if St is uniformly convergent and if St normalized by M 
is uniformly subconvergent on V^x,. 

How to check the uniform convergence of a scalar subdivision scheme is 
discussed in Subsection 6.2. The method is based on Theorem 13 from [4]. 
In the following we prove a sufficient condition for the uniform convergence 
of St normalized by M on V^, because to the best of our knowledge, such 
a theorem has not been published yet. Here we formulate and prove such a 
theorem for our specific case, namely for checking the smoothness of limits 
of a convergent scalar multivariate subdivision scheme, corresponding to a 
trigonometric polynomial in Z 1 . 

Theorem 17 Let tEZ 1 ,Mbeadxd dilation matrix. Then there exist 
trigonometric matrices T, Q of orders dxd and d 2 x d 2 respectively, satisfying 
(25) and (26). Moreover, ift(0) = m and 

||M* L || 00 ||^|| 00 < / u<l (45) 

for some positive integer L, then the subdivision scheme S T normalized by 
M uniformly converges on V^. 

Proof. The existence of T,Q, satisfying (25) and (26), follows from 
Theorem 10, as indicated in Subsection 6.1. 

To prove the second part of the claim, we will show that the sequence of 
vector-valued functions 

F k (x) = v(M k x - a)(M* k S T Vf) a , k = 1,2,..., 

with if as in the proof of Theorem 13', is a Cauchy sequence in (C(M. d )) d for 
any / G 4c- 
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Recall that ip is a continuous compactly supported function satisfying 
(30) and the interpolatory conditions (32), and that t, defined as in Lemma 
14, is in Z° and satisfies t(0) = m. 

By Lemma 14 with n = d, 

F k (x) = M* k J2 v(M k x-a)(S^f) a = M* k £ <p(M k+1 x-a)(S iId S$Vf) a . 



This yields that 

F k+1 (x) - F k (x) = M* k ¥(M k+1 x - a)((M*S T - S iId )S k S7f) a = 

M* k J2 ^{M k+1 x - a)(S f S k Vf) a , (46) 

where T = M*T — tl d - The entries of the matrix T are trigonometric poly- 
nomials in Z° because t E Z° as was mentioned above, and the entries of T 
are in Z° due to Theorem 10. Moreover, by Proposition 12, 

f (0) = M*T(0) - t(0)I d = M*mM*~ 1 - ml d = O d . 

Due to Lemma 3, it follows that for any g E 

S f g = SqS7 9 , (47) 

where Q is a d x d 2 trigonometric matrix. 
So, (46) may be rewritten as 

F k+1 (x) - F k (x) = M* k v(M k+1 x - a)(SQS7S k Vf) a , 

and 

\F k+1 (x) - F k (x)\ < CIIM^IUH^V^V/Iloo, (48) 

where C depends only on ip. 

It follows from (26) that VS^g = SqVg for any g E i^, and hence 

VS*V/ = S k V 2 f. (49) 
Thus we get from (48), in view of (45), 
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|F fc+1 (a;) - F k (x)\ < CIlM^IUH^IUII^vVlloo < 



Ci\\M* L \\J \\sk\\J |v 2 /l < c x |v 2 /|^. 



This yields that for all k, n e N 



n-l 



- F k (x)\ < d HvVlU^/iV < C 2 ||V 2 /||oo^, (50) 

i=o 

which implies that {Ffc} is a Cauchy sequence. Denote the limit vector- valued 
function by F which is, evidently continuous. Passing to the limit in (50) as 
n — > oo, we have 

(51) 



lim sup \F(x) - F k (x)\ = 0, 

fc^oo , 



Due to (51), 



lim sup 



F(x) - (M* k S*Vf)MM k x - a 



0, 



Substituting x = M k (3, (3 G Z d , we have 



lim sup 



F(M~ k p) - (M* k S k 7f)MP - a) 
The claim of the theorem follows now from (32). 



0, 
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